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Abstract 

The limit of small entropy production is reached in relaxing systems long after 
preparation, and in stationary driven systems in the limit of small driving power. 
Surprisingly, for extended systems this limit is not in general the Gibbs-Boltzmann 
distribution, or a small departure from it. Interesting cases in which it is not are 
glasses, phase-separation, and certain driven complex fluids. 

We describe a scenario with several coexisting temperatures acting on differ- 
ent timescales, and partial equilibrations at each time scale. This scenario entails 
strong modifications of the fluctuation-dissipation equalities and the existence of 
some unexpected reciprocity relations. Both predictions are open to experimental 
verification, particularly the latter. 

The construction is consistent in general, since it can be viewed as the break- 
ing of a symmetry down to a residual group. It does not assume the presence of 
quenched disorder. It can be — and to a certain extent has been — tested numeri- 
cally, while some experiments are on their way. There is furthermore the perspective 
that analytic arguments may be constructed to prove or disprove its generality. 



Key words: Nonequilibrium extensive systems. 
Pacs numbers: 75.40. Gb, 75.10.Nr, 02.50.-r, 05.20.-y 



1 Invited talk to be presented at STATPHYS 20, Paris. 

Paper dedicated to Heinz Horner on occasion of his 60 t/l birthday. 

2 leticia@physique.ens.fr 

3 jkurchan@enslapp.ens-lyon.fr 



Preprint submitted to Elsevier Preprint 



1 February 2008 



1 Slow dynamics 



1.1 Two complementary aspects: driven and aging systems. 



A universal feature of systems with slow dynamics is the extreme sensitivity 
of their time-dependencies to perturbations. Glassy systems which 'age' if left 
to relax (the correlations evolve slower and slower as time passes) may become 
stationary if an even small amount of power is pumped into them. Thus, for 
example, a 3D spin-glass having the characteristic aging autocorrelation curves 
(Fig. 3a)), becomes stationary upon applying a very weak, slowly evolving 
random field (Fig. 3c)). The reason for this sensitivity is the presence of 
flat directions in phase-space, indeed the very origin of the slowness of the 
dynamics. 

Remarkably, some properties remain unaltered even when the time-dependence 
has been dramatically changed by a small perturbation. Precisely these prop- 
erties are universal in slow systems, at least qualitatively. 

What this suggests is that we treat on an equal footing all problems with 
slow dynamics, whether they age or are stationary [1]. Their relevant common 
feature [4] (indeed, the definition of 'slow') is the smallness of the quantity 
that can be identified as an entropy production [5] in each case. Thus, we 
refer collectively as the 'small entropy production' (SEP) limit to the situation 
reached at long times and/or small external power input. 



1.2 Fluctuation- dissipation temperatures. 



In the SEP limit many systems are close to equilibrium. Such situation has 
been well studied in the past (see e.g. [5]). However, this is not the generic situ- 
ation, particularly in complex systems such as glasses, phase-separating fluids, 
moving objects that are just unpinning, etc. An example will at this point clar- 
ify things. Consider a system with variables 4>i(t), their correlation functions 



of 4>i to 

h=0 



Cij(t,t') = (£')), and the linear response Rij(t,t') = jf^jjj 

a kick applied to <pj. If we make a parametric plot of the integrated response 

t~\~t"W 

Xij(t + t w , t w ) = J dt 1 Rij(t + t w , t') (1) 



versus the corresponding correlation Cij(t + t w , t w ), the fluctuation-dissipation 
theorem is the statement that the plot gives a straight line of gradient —(3 = 



2 




0.2 0.4 0.6 0.8 1 

C(t + t w , t w ) 

Fig. 1. Integrated response versus correlation plot for a mean- field glass model. 

— 1/T. For an aging system we get instead, for large t w , a limit curve that looks 
qualitatively like Fig. 1. This figure has been obtained for the autocorrelation 
of the variables in a mean-field glass model [13], but similar plots have been 
obtained numerically for 3 and 4 dimensional spin-glasses [6], Lennard- Jones 
glasses with Montecarlo [7] and molecular dynamics [8]. Ordinary domain 
growth (as well as droplet models for spin-glasses) [9] also yield such a curve, 
with the peculiarity that the line to the left becomes horizontal. 

These systems age. As mentioned above, we can perturb them with small, 
nonconservative forces which render them stationary. Remarkably, at least in 
every model we know, if we now make the plots of Fig. 1 we get, even for 
small driving powers, the same dependence [2,3]. Despite the fact that the 
driven systems are stationary, they are well out of equilibrium. 

Such robustness of the plot in Fig. 1 suggests that the fluctuation-dissipation 
ratio might have a physical meaning. Let us define, for any two observables 
and for any two times, the effective temperatures as [10] 

uj\xeCij[uj, t w ) 

We have Fourier transformed the time-difference, as in experiments. In a sta- 
tionary system the t w dependence dissappears. 

1.3 Scales 

Equation (2) is just a definition. The scenario we shall discuss here consists of 
saying that for all pairs of observables and for each timescale there is a single 
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effective temperature. That is, 
%J f (LU,t w )-(3 e k { f (u,t w )~0 or (3tJ f (u,t w )~p e ff(u;,t w ), Vi,j(3) 

and 

/3 e// (u;i, * w ) - /5 e// (cu 2 , ~ if — = const. (4) 

in the SEP limit (t w — > oo and vanishing driving power). A system with a x 
vs. C plot like Fig. 1 (two straight lines) then has two temperatures (one of 
them being T). In systems with more than two temperatures the plot has a 
straight line of gradient — 1/T plus a piece of non-constant slope. 
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Fig. 2. Thermalisation of two subsystems. Integrated response vs. correlation plots 
for two uncoupled systems (dots) and for systems that are coupled (lines). 

1.4 Thermometry. Zeroth law. 

One can now ask if (2) really defines temperatures. To check this, one couples 
a 'thermometer' system simultaneously to the same observable of an ensemble 
of many independent copies of the system [10,11]. One can then show that 
the copies of the system act on the thermometer as a superposition of thermal 
baths with the corresponding effective temperatures and timescales. In order 
to measure a single temperature, one can construct a tuned thermometer (e.g. 
an oscillator) that responds to a single time-scale. 

To see the relation of FDT-temperatures with the zeroth law, one can couple 
the 'thermometer' as before, but alternatively to two observables. Then one 
shows that there is on average a net heat flow from the higher to the lower 
effective temperature. 
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1.5 Thermalisation of subsystems 

Within the present scenario, two systems with different effective temperatures 
if weakly coupled eventually 'thermalise' as follows: If the coupling is relatively 
strong [10], we have the situation of Fig. 2 (a direct solution of the equations, 
without assumptions): two (mean- field) glassy systems were defined so as to 
have different effective temperatures when uncoupled (dotted lines). If they 
are made to evolve coupled, the effective temperatures equalise (full lines). 

If the coupling is weak [11], the system preserves essentially the tempera- 
tures of its constituents, but eventually rearranges them in widely separated 
timescales. Thus, the combined system has three temperatures, the bath tem- 
perature plus other two, each acting on a separate timescale. 

1.6 Auxiliary thermal baths 

An interesting tool for the study of a system with more than one 'natural' 
temperature is to couple it to a slow, weak auxiliary bath [1 1] of temperature 
T*, which, as we shall see, plays the role of 'a field conjugate to the natural 
temperatures' [12]. The outcome, within this scenario, is that the auxiliary 
bath fixes the timescale at which its temperature happens in the system. In 
particular, an aging system with a single T e ^ > T will become stationary 
if T* > T e ff, and be hardly affected if T* < T e ^ . More generally, an aging 
system with multiple effective temperatures should become partially station- 
ary (for all the timescales with T* > T e ^), but would still have aging for 
timescales with T* < T e ^ . Simulations for the three dimensional Edwards- 
Anderson model (a system which seems to have many temperatures [6]) are 
in this sense very encouraging (see Fig. 3). 

1.1 Reciprocity 

A surprising feature of the present scenario is that in the SEP limit there are 
reciprocity relations (A(t)B(t')) = (B (t) A (t')}, and similarly for the responses. 

It is remarkable that reciprocity relations hold in a situation in which FDT is 
strongly violated. This is all the more surprising in an aging case in which the 
system is not even stationary. The interest of these relations is that they are 
relatively easy to measure in a simulation or in an experiment. 

Figure 4 shows a numerical check for two 3D spin glasses. C\\ and C22 denote 
the autocorrelations of systems 1 and 2 (which have coupling independent 
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Fig. 3. Autocorrelations of 3D Edwards-Anderson model. From top to bottom: a) 
pure relaxtion, b) same but coupled to a weak slow bath with T* = 10, c) same but 
coupled to a weak slowly evolving random field (a weak slow bath with T* = oo). 
In b) aging has partially dissappeared, while in c) it has dissappeared completely. 
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constants, on average half as strong in one system than in the other). Despite 
the assymmetry between the two systems (cfr. the left figure), the mutual 
correlations Cyi and C21 tend to equalise more and more for longer waiting 
times (right figure). 



1.8 Emergence of macroscopic temperatures. 



A very encouraging feature of the effective temperatures, such as they come 
out in the solution of some analytic models [13], is that they remain non-zero 
in the limit of zero bath temperature. In other words, the temperatures stay 
non-zero even when Boltmann's constant is negligible. This is most welcome: 
if these effective temperatures are to be relevant for macroscopic cases such 
as colloids, sheared foams and (perhaps) granular media, they must be huge 
in units of the Boltmann constant, and be related to macroscopic structures. 
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4. Parametric plots of Cxi vs. C22 (left) and C12 vs. C21 (right) for t w = 
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Fig. 

(diamonds), t w = 10000 (crosses), t w = 30000 (squares). Note that the points in the 
right graph approach the diagonal for larger t w . 



2 Analytical results. Symmetry breaking 



Let us outline very briefly the analytical arguments showing that this scenario 
'closes'. We assume the system has Langevin dynamics: 

dV{6) 

- m^i - Mij^j — — = r o 0j - rji (5) 

d(pi 

where rji are uncorrelated Gaussian white noises with variance 2r T and r 
is the strength of the coupling to the 'white' bath. One can also add weak 
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nonconservative forces. The masses m 8 can in particular be zero. 

We can encode the correlations and (causal) responses in the 'superspace' 
order parameter (see e.g. [14]): 

Qij(l,2) = Cij{t u h) + (0 2 - h) [02 R^itiM) - di Rjifah)} . (6) 

where a , a are Grassmann variables, and we denote the full set of coordinates 
in a compact form as 1 = dl = dtd0d0, etc. 

One can obtain a general perturbative expression for the correlations and 
responses ([11], see also [14]): 

m t ^ + T D^ g ij (l,2) + M ifc Q fcj -(l,2) = <J(l-2)(J - + [E(8)Q] ij (l,2) (7) 

where <g> stands for convolution and matrix product. The function S is ob- 
tained from all diagrams with propagator Q with two (amputated) legs that 
cannot be disconnected by cutting one more line. Note that we have nowhere 
assumed there is quenched disorder, and that this expansion can be done for 
any model, independently of the dimensionality. 

We now make a separation between 'fast' Q F and 'slow' Q evolution Q = Q F + 
Q. In the 'fast' regime of finite time-differences, correlations and responses can 
be shown [4] to satisfy FDT, and are time-translation invariant. In the SEP 
limit one can obtain a separate equation for the 'slow' part: 

M ik Q kj (l, 2) - %j T %j £[QU1, 2) - hj T kj Q ik (l, 2) = 'small' (8) 

Here 'small' stands for terms that vanish in the SEP limit (see below), and 
S[Q]ifc(l,2) is a functional of Q that is also obtained diagrammatically. The 
Qkji Qkfi dfj, dfj A are parameters to be determined from (8) and the equations 
for Q. 



2. 1 Symmetries 

In general, to the extent that one is allowed to neglect the 'small' terms in the 
SEP limit, Eqs. (7) are invariant with respect to any change of 'coordinates' t a , 
a , and a (a = 1, 2, ...) with unit superjacobian [15]. This is a large symmetry 
group, including in particular the following: 
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i. Ordinary time-reparametrizations. 

t a ^K(t a ), 9 a ^K(t a )9 a , 9 a ^9 a . (9) 

These reparametrizations were already considered by Sompolinsky and Zip- 
pelius in their seminal work [16]. Under this subgroup the order parameter 
Q preserves the form (6). 

ii. Time-reversal. 

t a —> —t a — (3*99 , 9 a — > 9 a , 9 a — > 9 a . (10) 

These reparametrizations only preserve the form (6) if they act on a Q such 
that Rij(t - t') = f3*d t >C i:j (t - t') G(t - f). 
iii. Boson-Fermion symmetry 

A group of symmetry transformations which, unlike the previous two, 
exchanges bosons and fermions is 

K(t a )^K(t a ) + K , 6 a ^9 a + e, e a ^e a + i, + ^—9 a (11) 

K (t) 

for any K(t), K , e, e, e (the latter three, odd-Grassmann parameters). 
2.2 Temperature fixing as symmetry breaking 

Let us now briefly show how this symmetry considerations prove the consis- 
tency of the ansatz. We shall for brevity only mention here the case of a single 
effective temperature T e ^ for the slow dynamics. 

Equation (8) (neglecting its r.h.s.) will always admit a solution invariant with 
respect to some residual subgroup of the group of transformations with unit 
superjacobian. To get a non-constant solution, we must then break (i), by 
choosing a time-reparametrisation K(t). One can easily check that fixing an 
effective temperature now corresponds to breaking (iii) to a three-parameter 
subgroup with fixed K(t) and with the restriction e = T e ^e in (11). One can 
further keep 'time reversal' (10), but as applied on the reparametrised time 
K(t), and with T* = T e ^ . 

This may seem rather strange, but in fact it is the abstract way of describing 
the usual temperature fixing by a thermal bath. In any purely Hamiltonian 
system (T = 0), Eqs. (7) are invariant under the group (11) with K(t) = t 
and four parameters Kq, e, e and e. If we now switch on a thermal bath of 
temperature T (we make Tq > 0), the group (11) is broken down by the 'bath' 
term proportional to T to the residual three parameter subgroup [17] obtained 
by restricting e = Te. 
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2.3 Sensitivity and the matching problem 

The neglect of the 'small' terms in the right hand side of equation (8) brings 
about the invariance with respect to super-reparametrisations. However, the 
true equations do not have so many symmetries. In order to see what is the 
situation, we note that equation (8) is of the form: 

[ (Super) Reparametrization Invariant ] = 'small', non-invariant (12) 



If we find one solution to this equation by setting the l.h.s. to zero, any 
super-reparametrization is also a solution. If the solution to l.h.s.=0 breaks 
reparametrization invariance, one obtains thus an infinite number of solutions. 
In order to fix the correct one, one must take into account the 'small' terms in 
the r.h.s. This requires a calculation going beyond the zero entropy production 
limit (small but finite forcing or long but finite waiting times), and this is a 
much harder task, indeed a problem that remains to be solved. 

On the other hand, apart from being a technical nuisance, this has a physical 
meaning. Because terms that are arbitrarily small decide which reparametriza- 
tion is the good one, systems having non-trivial solutions to (8) are sensitive 
to perturbations that are arbitrarily weak, a general feature of systems with 
slow dynamics. 



3 Perspectives 

Given the number of not entirely obvious relations this scenario proposes, 
the least one can now say is that it has heuristic value, even if the fact of 
having been inspired by mean-field spin-glass theory is sometimes viewed as 
an original sin. There has been recently quite a lot of numerical activity to 
test the predictions, with encouraging results. The experimental exploration 
of structural glasses from this point of view is also now beginning. 

The solution is consistent in general, since it corresponds to the breaking of 
a symmetry that is a general property of slow evolution down to a residual 
subgroup (having the physical meaning of partial thermalisations). Just as 
in any partial symmetry breaking situation, one cannot with such arguments 
alone assure that in some model this symmetry will not be broken further, or in 
a different way. One may perhaps hope that real proofs concerning the generic 
situation in the SEP limit of extended systems will be one day available. 
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